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Recently, the stochastic series expansion (SSE) has been proposed as a powerful MC-method, 
which allows simulations at low T for quantum-spin systems. We show that the SSE allows to 
compute the magnetic conductance for various one-dimensional spin systems without further ap- 
proximations. We consider various modifications of the anisotropic Heisenberg chain. We recover the 
Kane-Fisher scaling for one impurity in a Luttinger-liquid and study the influence of non-interacting 
leads for the conductance of an interacting system. 
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I. INTRODUCTION 

In general, integrable one-dimensional models show 
an ideally conducting behavior in contrast to most real 
three-dimensional materials.^ There has been intensive 
investigation of the influence of the integrability on the 
conductivity for various model systemsi^ii^ The con- 
ductance of conducting, nearly one-dimensional devices 
is, on the other hand, of substantial experimental in- 
terest. Over the last years it has become possible to 
fabricate mesoscopic devicesf^ such as carbon nanotubes 
which can be viewed as a realization of systems with bal- 
listic transport properties. Therefore, the computation of 
dynamical transport quantities has received considerable 
interest. 

A basic approach for the study of the con- 
ductance has been in the past via the bosoniza- 
tion of appropriate model systcmSfSi^ valid in the 
low-temperature limit. Numerical studies have so 
far involved the density-matrix-renormalization-group 
(DMRG) techniqueSiia and Monte Carlo (MC)iiii^ In 
the case of Ref. |^ a reduced set of states was used in or- 
der to evaluate the dynamics, in Ref. llOl a phenomeno- 
logical formula by SushkoA*i^ was used to compute the 
conductance (see also Ref.0). The simulations by Refs. 
ITllIT^ use an effective bosonized Hamiltonian as a starting 
point. 

Here we will discuss how to obtain the conductance 
with quantum-Monte-Carlo (QMC) on the original lat- 
tice Hamiltonian. For this purpose the conductance will 
be calculated on the imaginary frequency axis. We will 
show that a reliable extrapolation to zero frequency can 
be performed at finite but low temperatures. We will 
thus obtain approximation-free results for the dynamics 
of inhomogeneous quantum-spin systems at low but finite 
temperatures, within a well defined numerical accuracy 
defined by the statistics of the MC-sampling and the ac- 
curacy of the zero-frequency extrapolation. 

By the Jordan- Wigner transform a one-dimensional 
spinless fermionic system can be mapped to a hard-core 
boson model. Hence, it is possible to calculate the con- 
ductance for a fermionic system in a bosonic one. This is 



vitally important since boson models can be easily ana- 
lyzed by Monte Carlo simulations where the sign problem 
is absent. However, for an evaluation of the conductance 
one requires a highly efficient simulation method which 
performs well at low temperatures. RecentlypSi such a 
powerful method has been proposed: the Stochastic Se- 
ries Expansion (SSE). In this paper we will compute the 
conductance in a hard-core boson lattice model by the 
aid of this new method. 



II. EXPLANATION OF THE METHOD 
A. Definition of the Conductance 

We consider the anisotropic xxz-Hamiltonian 

N-l 
n=l 

where the = 5*^^ ± iS^ are the raising/lowering op- 
erators for spin-1/2 Heisenberg spins. The spin-current 
operator j„ at a given site n follows from the continuity 
equation and is given by (see e.g. Ref. 0) 

jn = iJxG {S^Sj^^i — S'„ S^^i) /{2h). 

As a perturbation we will use a local "voltage drop" . In 
the hard-core boson notation this corresponds to a step 
in chemical potential — at site m — which is equivalent to 

in terms of the Heisenberg-spin operators. 

The conductance is then defined as the dynamical re- 
sponse of the current operator at site x to the voltage 
drop at site y: 

g:=lim Re-/ e^\[jx{t), Py]) dt. (1) 

For open boundary conditions (OBC) the relation 
i[H,Px\ — hjx holds and a partial integration of 
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yields: 



Re 



Using Re(a6) = ReaRe5 — Imalmfo in the above equation 
gives two contributions to the conductance. With the 
definition of the generalized Drude weight for two oper- 
ators A and B : 
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where AA ^ A — (A) one can show using the Lehmann 
representation^^ that the first contribution (ReaRe6) 
reads 



{{jxjv))'R.e{-iz) = {{j^jy))7T6{Rez). 



(2) 



The second factor in the expression ([SJ [namely, 
Re(— gives rise to a delta- function, such that (0) 
should be the dominating contribution to g. The first fac- 
tor of expression Q (namely, {{jxjy))) is closely related 
to the Drude Peak D = {{JJ))/N where J = E„J" is 
the total current operator. 

From the discussion of the Drude Peakii we know 
that under OBC's the Drude Peak is zero [because D 
can be written as the response to a static twist which 
can be removed by a gauge transformation of the form 
exp(z^„n5'^) (see Ref. [TJ) ] whereas it is non-zero 
for periodic boundary conditions (PBC's). In our case 
{{jxjy)) is zero under OBC's [use a gauge transform 
exp(i^^^^ S^)]. Under PBC's we find, because of trans- 
lational invariance and because of the continuity equa- 
tion, that {{jxjy)) does depend neither on x nor y. This 
implies D — N{(jxjy)). Since the Drude peak is finite 
(at least for the models we are interested in) we con- 
clude that expression ||2Jl vanishes even under PBC's in 
the Thermodynamic limit. So we obtain 



g = Re[{zh)-^]Re / e'^' {[jx{t), jy])dt. 
Jo 



(3) 



Restarting from Eq. we may — again by partial 
integration — arrive at another formula. 
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[Px{t),Py])dt} 



The first term in the square brackets does 
contribute — as the potentials Px and Py commute— 
if we restrict ourselves to Rez = we obtain: 
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e''*{[Px{t),Py])dt 



not 
-and 



(4) 



The latter formula is especially useful for MC-simulations 
as it allows the computation of the conductance in terms 
of the diagonal 5'^-S'^-correlators (under OBC's). 
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FIG. 1: The conductance of the xy chain at T = 0.01 versus 
frequency on the imaginary axis for various distances {x — y = 
2n, n G N) between voltage drop (y) and current measurement 

(x). 



According to its definition as it is given by Eq. 
the conductance might in principle depend on the actual 
choice of the positions of the voltage drop y and the cur- 
rent measurement x. Here, we point out that in the limit 
z —^ this is not the case. In a rather general situation 
one can show (using the continuity equation) that the 
right hand side of Eq. |^ gives the same result for any 
choice of x and y. (see Appendix El) 

It is instructive to consider the free fermion case en 
detail. We denote the — formal — dependence on x and y 
by corresponding subscripts. Of course, in a translational 
invariant system g^y depends only on the difference x — y. 
The conductance g as a function of w = Imz (here and in 
the sequel Rez = 0) is plotted for the free fermion case 
in Fig. n] One sees that the conductance in the limit 
Lu approaches the universal value e^/h. Here we 
emphasize that a spatial separation of voltage drop and 
current measurement leads to an exponential decrease in 
g{uj) at small uj. Therefore, we will restrict our attention 
to the cases \x — y\ < 1 for the rest of the paper. 



B. Technical details of the MC-method 

We now turn to some technical details of our simu- 
lations. The second formula for the conductance, Eq. 

lends itself to a study with Monte Carlo simula- 
tions (it requires OBC's). At the Matsubara frequencies 
ujm = 2TrM{Ph)~^ , M € N we may use the equivalent 
expression 



g{ujM) = -ujAi/hRe / {PxPy{iT))e 



'dr. 



We employ a standard QMC method (SSE) to com- 
pute the conductance. Since Px is diagonal in the S^- 
Basis, the simulation of {H^PxH^~^Py) can be easily 
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performed with the help of the SSE.^^-^^'^^ Here, L is the 
approximation order. One may simply obtain {P^Py{iT)) 
as a linear combination of the terms {H'^PxH^^^Py) 
with binomial weight factors B{t^ k). We found it conve- 
nient to assume a Gaussian distribution for the B{t, k) 
instead of a binomial one, because the former is easier 
to evaluate. The error that we introduce by this replace- 
ment is smaller than the statistical error if L > 100. 
(Note that in our simulations L is typically of the order 
of 10'' - 10^.) 

What remains to be done in order to get 

g{0JM) = -^Al/fl I C0S{UJ MT){PxPy{iT))dT 

Jo 

is an integration in the final step. We performed it with 
the Simpson rule and a grid of 800 r-values. 

We are now left with the standard problem of extrapo- 
lating g{uj) from the Matsubara frequencies com to uj — 0. 
Unfortunately, the spacing of the Matsubara frequencies 
is linear in T, so our method becomes unstable when we 
increase T. 

To see that an application of our MC-method makes 
only sense at low temperatures we compare it to a sim- 
pler method: exact diagonalization. Fig. |21 shows g{uj) 
for various system sizes at T = Jx/^b- One sees that 
convergence with N is rather fast (at high temperatures). 
Hence, one can determine g^uj) for uj > Jx/fi with exact 
diagonalization. To compute g{Lo) at some lo < Jx/fi 
with MC-methods one needs to work at a temperature 
T < Jx/ {'i-T^ks), and even then exact diagonalization is 
preferable as it yields g{oj) in a continuous interval rather 
than on a discrete set of points. Hence, at high tem- 
peratures the MC-method is inferior to a simple exact 
diagonalization. 

In our simulation we make one "MC-sweep" between 
two measurements which consists of one diagonal update 
and several loop-updates-'^^ between two measurements. 
We are able to run approximately 10^ sweeps. 



C. Test with Jordan- Wigner 

If J2 = then g{u}) can be exactly evaluated — for ar- 
bitrary system size — not only at w = (see below). 

We can exploit this fact in two regards: Firstly, we test 
our MC-method by comparing it with the exact curve 
obtained by Jordan- Wigner. The result can be seen in 
Fig.El 

Secondly, we can test which frequencies (and hence 
which temperatures) we need such that a linear extrap- 
olation can be carried out without introducing a larger 
error than the statistical one. One sees also which system 
sizes are needed to determine g{ijj) without detectable fi- 
nite size error. Our conclusion is that our method works 
for T < 0.02. At T « 0.01 a system size of N k 200 - 300 
is appropriate. 



III. THE CONDUCTANCE IN VARIOUS 
SYSTEMS 

A. Results for a dimerized Jordan- Wigner-chain 

One may derive a simple analytical result for the con- 
ductance in the free fermion case. Here, we consider the 
slightly more difficult case — but also more interesting as 
the system has a gap^S — of a dimerized chain with mag- 
netic field B, i.e., = and the hopping parameter 
alternates: {Jx)2n,2n+i = Ji and ( J2;)2n+i,2n+2 = ^2- 
The energy dispersion is given by 

Ef =B±^Jl + Jl + 2 Ji J2 cos(2fc) / 2 

and we assume a positive dispersion Ek = E^ if fc G 
[0,7r/2] U [37r/2,27r] and Ej, = E' else. The gap Eg is 
Eg = E^^^—E~^^ at fc = ±7r/2. The evaluation of the for- 
mula for the conductance Eq. ^ is in principle straight- 
forward, (see Appendix IbI) One obtains the compact 
result: 

g = 1^ [tanh(So/3/2) - tanh(^,/3/2) 

- tanh(ii;+ 2/3/2) + tanh(£;;/2/3/2)] . (5) 

In the case of zero magnetic field this reduces to: 

g = [tanh(^o/3/2) - tanh(S+ ^/S/S)] . 

The T — 0-value of the conductance as a function of 
magnetic field is quantized: It is zero if \B\ is smaller 
than the zero field gap E^^^ or larger than the zero field 
band width Eq, and it is 1 between these values, and 
precisely at these values it is 1/2 (all values in units of 
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FIG. 2: Conductance g as a function of tj. The arrow indicates 
the result of the Thermodynamic limit. 
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FIG. 3: Monte Carlo data (symbols) in comparison with the 
exact Jordan- Wigner result (dotted lines). The curves for 
T = 0.02 are offset by 0.1 for clarity. (We use OBC's.) 



FIG. 4: Monte Carlo data (symbols) for various interaction 
strengths at T = O.OlJ^/kB (using OBC's and 2 • lO'^ Monte- 
Carlo sweeps). Shown is the conductance as a function of 
imaginary frequency. 



B. Comparison with the Apel-Rice- formula 

At low temperatures the xxz-chain can be described 
by a Luttinger liquid. For this model the conductance 
was first obtained by Apel and Rice in the eighties>2i2i 



<7ApclR,icc 



h 2{Tr-e)'' 



(6) 



where cos9 — Jz/ Jx- This formula may be derived from 
Eq. |@J if we use concrete expressions for {S^S^{iT)) 
which are available from conformal field theory^^. Fig. 
0]shows our QMC-results for g{uj) on the imaginary axis 
for the xxz-model, in Fig.|21we display a comparison, as 
a function of Jz, between the g{uj = Q) extrapolated from 
Fig. 21 and the exact Bosonization result, Eq. ©. 

We note, that the statistical error of the QMC-results 
presented in Fig. does not increase much with the 
parameter Jz- This is due to our using the "directed 
loops" as described in Ref. [^J The choice for transi- 
tion probabilities which was proposed there makes the 
SSE-algorithm more effective. Here the improvement is 
remarkable. 



C. System with one impurity 



Now, we consider the Hamiltonian 

i.e., we add an impurity in the middle of the system. This 
model was first studied in a paper by Kane and Fisher 
Later, this kind of model received considerable atten- 
tion by other authorsiiSiS^iSiSSiSi By an RG approach 
Kane and Fisher found that the perturbation Bimp is 
relevant for repulsive interactions (i.e., Jz > 0) . This 
means that at zero T the chemical potential anomaly 
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FIG. 5: Conductance of the a;a;2:-spin chain at T = 0.01Ja;/fc_B 
in comparison with the exact result by Apel and Rice for 
a Luttinger Liquid. The Monte-Carlo data are obtained by 
extrapolating the results of Fig.|l]to u — 0. 



"cuts" the system into two halves, such that the conduc- 
tance is zero. This result cannot be directly confirmed by 
Monte Carlo methods since these are necessarily finite- 
temperature methods. 

Fortunately, the scaling behavior (with temperature) 
of the conductance is also known. For K = 1/2 one 
may derive an exact formula for the conductance by a 
refermionization techniquc:^'^"' 



r>2 



1/2- 



(7) 



where "0 is the Digamma function. So we can compare 
our MC-data once again with an exact result. In Fig.lHI 
we present two different QMC-results for the conductance 
on the imaginary axis of the Heisenberg-chain with one 
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FIG. 6: Conductance of the Heisenberg-chain (A'^ — 400 sites) 
with one impurity at T = O.OlJx/fcs for different Impurity 
strengths and positions of voltage drop and current measure- 
ment (using OBC's). A possible extrapolation to = is 
proposed by the dashed lines. 



impurity, for different impurity strengths. 

For the upper set of curves in Fig. the position of 
voltage drop and the position of the current measurement 
are x = y = N/2; and for the lower set of curves they 
are x = N/2, y = N/2 — 1. The curves are not as 
smooth as the ones in Fig. 01 so we use a quadratic fit 
from the first three Matsubara frequencies instead of a 
linear extrapolation to estimate g{uj = 0). We also note 
that the curves with x — y — 1 are better suited for 
extrapolation than those with x = y because the slope 
at = is smaller. The statistical error is less than 
one percent. The result from the extrapolation is given 
in Fig. [7| along with the exact curves from Eq. Q . We 
used system sizes of = 400 for T > O.OlJx/kB and 
iV = 800 for T = 0.005 J^/fcij. We performed 2 • 10^ MC- 
sweeps. The error bars are smaller than the symbol size, 
so the error that we see in the figure is mainly due to our 
extrapolation method [and to possible logarithmic finite- 
temperature-corrcctions to Eq. Q]. We see that the 
quadratic fit tends to underestimate the correct value. 



D. Inhomogeneous systems 

Several years after the publication of the Eq. ® by 
Apel and Rice it was generally agreed uponSSiSS that it 
does not refiect the (correct) physical behavior one would 
encounter in an experimental realization. Experiments 
are never performed on a closed system but on one cou- 
pled to reservoirs which make it possible for the particles 
to leave and enter the system. These reservoirs can be 
modeled by attaching two leads consisting of infinite non- 
interacting spin half-chains to our model. 
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FIG. 7: Conductance of the Heisenberg-chain with one impu- 
rity at various temperatures (A'^ — 800 for T = O.OOSJx/fcs 
and A'^ = 400). MC-data (symbols) are drawn in comparison 
with the exact formula (solid lines) by Weiss et ali— for a 
Luttinger Liquid at the isotropic point K = ^ (i.e., Jz = Jx)- 

The complete Hamiltonian reads then: 

(JV+Wj-3)/2 

H = HxxziJz = 0) + ^ S^S^^i, (8) 

n={N-Ni + l)/2 

i.e., the interaction is confined to a small region in the 
middle consisting of Nj sites. This appro ach has been 
followed by many authors, e.g. Refs. ITollsollsllls^ . Gener- 
ally, the presence of leads yields a conductance which is 
independent of Jz jSSiSS namely, 

9 = e'/h 

in sharp contrast to Eq. ((BJ. The non-interacting semi- 
chains — which we call leads — play the role of reservoirs. 

We note that the behavior of g{uj) depends on the par- 
ity of Nj a fact which was already reported in Ref. 10. (A 
similar effect was also found in the Hubbard model. '^■^) In 
the following we will only consider the case of Nj odd. 
A detailed discussion of Nj even/odd and a comparison 
with Ref. 01 will be presented elsewhere. 

The "natural" choice for current measurement and 
voltage drop would be at the two ends of the interacting 
region. But here is caution advised. From Fig.^ — which 
is again for the free Fermion case — one learns two things: 
Firstly, g does not depend on the choice of x and y. Sec- 
ondly, if X and y are some distance apart, convergence 
with Lu becomes slow, hence one needs to go to lower T 
and larger system sizes if one wants to extract g{uj = 0) 
reliably. A simple phenomenological explanation for this 
is the following: If the place of the measurement is far 
from the voltage drop the particles have to travel a long 
distance and hence one has to wait a long time, before 
one can determine g. We can now place both the voltage 
drop and our current measurement at the middle of the 
system, but this should not help much. The problem is 
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FIG. 8: Conductance, in units of e^/h, of the Heisenberg 
chain with non- interacting leads, see Eq. Q. The system size 
is = 320. The values for temperature T and size of the 
interacting region Ni are given. (We use OBC's, 10^ MC- 
sweeps, and x = y for the voltage drop.) 



that the particles stiU have to travel a long distance un- 
til they see the leads. Hence g{oj) will be unaffected by 
the introduction of leads if cu is sufficiently large. This 
expectation is confirmed by the QMC-data presented in 
Fig. |S1 In this figure the exponential decay of g{u!) at 
small CO is apparent, and the curves illustrate clearly that 
this decay is induced by the length scale Nj — because 
it becomes stronger with increasing Nj. It is this decay 
that prevents us from discussing larger Nj. If we increase 
Nj the decay becomes stronger, hence we need to eval- 
uate g{Lo) for more (and smaller) frequencies in order to 
extract g{Lu — 0) reliably. But smaller frequencies are 
only available at smaller temperatures. As we cannot 
decrease T much below O.OlJx/fcs we restrict ourselves 
to Nj < 20. If we considered a system with Nj — 200 
(at T = O.OlJx/kB) we would not see any difference from 
a system without leads, because the difference occurs at 
small Lo. 

The data presented in Fig. |S1 clearly indicates an up- 
turn of g{Lo) for LO ^ 0, indicating that g{u! — 0) 
is unaffected by the interaction in the low-temperature 
limit. Our result may, however, not be totally convinc- 
ing, since we can only analyze relatively small interacting 
regions. One might argue that the enhanced conduc- 
tance is not due to the leads but simply to the reduced 
"mean" interaction — which is close to zero as only few 
sites interact. To invalidate this argument we considered 
another model. We have performed QMC-simulation of 
a system where we attach a lead only at one side such 
that we obtain a chain which is non-interacting in one 
half and interacting in the other. For this system we 
found no deviation at all imaginary frequencies from the 
situation where the interacting region extends over the 
whole chain, even though there are as many interacting 
as interaction-free bonds. 



E. Spin-Hamiltonian with third-nearest-neighbor 
interaction 

Monte Carlo simulations allow the inclusion of long- 
range hopping, and thus breaking the integrability of the 
pure xxz-Hamiltonian, as long as the resulting system is 
not frustrated. We thus consider a Hamiltonian with a 
third- nearest-neighbor interaction: 

H = Hxxz + Jx3 [ + S„ S^^^)/2 

n 

+ J,/J,5,:5,:+3]. (9) 

For simplicity we assumed that the anisotropy is indepen- 
dent of the hopping range (i.e., Jzs = JxsJz/ Jx)- In this 
context we emphasize, that the long-range hopping in the 
spin system does not transform under Jordan- Wigner to 
a long-range hopping in a fermionic system but to a more 
complicated four-sites operator. 

In general, adding a new term to the Hamiltonian 
changes the current operator which is defined via the 
continuity-equation, Vj := (j„+i - j„) = i[S^_^,i, H]/h. 
In a one-dimensional system with OBC's the continuity- 
equation is solved however by the relation jx = i[H, Px]/}i 
(see Sec. Ill A|l such that Eq. Q still applies. Nonethe- 
less, it is useful to look at the current operator for this 
case. It reads: 

jn = jn,l + in, 3 + jn-1,3 + jn-2,3, 

where jn,k = iJxB {S:^S~^^. - S'S^j^^) /{2h). If we com- 
pare it with the current operator of the xxz chain we 
see that the long range hopping Jx^ gives rise to three 
additional terms which are analogous to the first one. 

We compute the conductance as a function of the hop- 
ping amplitudes Jxz/Jx and present the result in Fig. 
|51 If Jx3 — the conductance is of course given by 
the Apel-Rice- result Eq. ©. However, if Jx^ » Jx 
we may eventually neglect the nearest-neighbor-hopping- 
term such that we end up with three uncoupled chains. 
Thus, we conclude that the conductance will grow to- 
wards three times the Apel-Rice-result when we increase 
Jx3- ^From the figure we see that the crossover between 
these two values is shifted to smaller values of Jx3 when 
the anisotropy is increased. (In fact, at the isotropic 
point the increasing of g is barely visible.) 

In conclusion we have developed a QMC-technique 
which allows the evaluation of the DC-conductance for a 
wide range of non- frustrating quantum-spin chains at low 
but finite temperatures. We have presented several strin- 
gent tests for this technique, like the Kane-Fisher scaling 
for the conductance through an impurity in a Luttinger- 
liquid. 
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FIG. 9: Conductance of the Hamiltonian Eq. Q (with 400 
sites, T = 0.01) as a function of Jx3 for various anisotropies. 
squares: quadratic fit from the first three Matsubara frequen- 
cies; diamonds: linear fit from the first six frequencies. (We 
use OBC's, 10^ MC-sweeps, and x = y + 1 for the voltage 
drop.) 



APPENDIX A: PROOF THAT THE RIGHT 
HAND SIDE OF EQ. ^ IS INDEPENDENT OF x 
AND y 

Here we provide a general argument which relies on 
the (physical) assumption that some linear response func- 
tions are finite in the Thermodynamic limit. 

Theorem: In a Spin system the conductance gxy = 
g does not depend on x and y if the linear response of 
an operator to a perturbation Pm is hounded in the 
Thermodynamic limit. 

Proof: 

First we consider lim^^o z e*^*([S'^(t), Pm])dt. This 
expression corresponds to a plateau value of the response 
function (p{t) := i{[Pm, S^{t)]), i.e., 



lim z 



{[S:,{t),Pr.])dt= lim ^{t). 

t — >oo 



The response function may be written in the following 
way (by Kubo's identity)^ 

(p(t) =/3(F„„zLe*^*5;j ^(t). 

Here (•,•) is the Mori scalar product [for operators 
A and B: (i{B,A) = /^^ TrSt exp(-Ti?)Aexp((T - 
/3)iJ)(iT/Trexp(— /3i/)] and L is the Liouville operator. 
To prove that Lp{t) Q as t oo, it is sufficient to show 
that $ as a function of t is bounded. 

But this is just the assumption that we made in the 
statement of the theorem because {S^{t),Pm) represents 
the linear response of the operator to the perturbation 
p 

Finally, we can prove our main assertion. We want 
to show g{x,y) = g{x',y') \/x,y,x' ,y' . This follows 



easily from g{x,y) = g{x + l,y) Vx,y and g{x,y) = 
g(x,y -\- 1) yx,y. We will only consider only the sec- 
ond equality [the proof of the first equality is analogous 
by the symmetric structure of Eq. Q]. 

Using Eq. Q and our previous result we see: 



^{x,y)^g{x,y+l) = \imlmz / e^^*([5^\i(t), P,])dt = 



Q.E.D. 



APPENDIX B: DERIVATION OF EQ. dill 

For eigenvalues En and Ek the respective one-particle- 
eigenstates will be denoted by \n) and the annihila- 
tion operators by c„ and Ck, and the occupation numbers 
by n„ and rifc. 

For the current operator we find: N{k,jxk) — evk '— 
e/h^ and {-k,j,k)=0. 

In a free Fermion system one can derive a simple ex- 
pression for {A{iujM)B) = dTe''^'"'^{AB{iT)) [with 
iojM — Lu + iS = 2TTiM/ (HP)] when A and B are one- 
particle operators (i.e., A = J2n k ^nkclck): 



ihwM + En - Ek ' 



where f, - sinhif3{E„-E^)/2) _ , , .^ 

Wtiere Jkn - 2cosh(l3E„/2)cosh((3Ek/2) " \\^k) ~ {Jin)} ■ 

{A{z = Lj -\- iS)B) may then be obtained by analytic con- 
tinuation. In our case A and B are local current opera- 
tors; for the conductance we use Eq. 



Ij {n,jxk){kjyn) 

I™ TT — ~ 



uj ^ ihwM + En - Ek 

n^k 

--^{^Jxk){k,jyn) 



{hw + En - EkY + 5^ 



fkn- 



For (5 



{hw + En - Ekf + 5'' 



TrS{hLU + En - Ek). 



In the continuum limit we replace N ^ ^ / ~ 

f dEk 

Performing the integration over the variable n and then 
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taking a; ^ yields Note that the conductance does not depend on the 

energy dispersion but only on the band width Eq and 



2hJ 2tt^ ^ '^2cosh^(f3Ek/2) 



2cos\i\fiEkl2) 
dk ftwfcSgn(fc)/3/2 f^='' dE(3/2 



h J_,2TTAcosh\(3Ek/2) h Jk=o 2cosh\f3Et,l2) 

g2 .So/3/2 1 g2 ^ 
= — / dx 2 1" oT / 2 

•^^^/2/5/2 '^osh (x) 2/i Je„i3/2 cosh (a;) 
= |^[tanh(£;o/3/2)-tanh(i;^/3/2) 



- tanh(£;+ 2/?/2) + tanh(£;;,2/?/2)] . 
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